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FUEL-OPTIMAL TRAJECTORIES FOR AERO ASSISTED 
COPLANAR ORBITAL TRANSFER PROBLEM 
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ABSTRACT: The optimal control problem arising in coplanar, orbital transfer 
employing aeroassist technology is addressed. The maneuver involves the 
transfer from high Earth orbit to low Earth orbit. A performance index is 
chosen to minimize the fuel consumption for the transfer. Simulations are 
carried out for establishing a corridor of entry conditions which are suitable 
for flying the spacecraft through the atmosphere. A highlight of the paper is 
the application of sin efficient multiple shooting method for taming the 
notorious non-linear, two-point, boundary value problem resulting from 
optimization procedure. 


NOMENCLATURE 

A : constant = Sp /2m 

S 7 

A : constant = C Sp H /2m 

1 dO s a' 

A constant = C Sp H /2m 

2 LR b a' 

a : ratio of R /R 

d da 

a : ratio of R /R 

C c a 

b : constant = R /H 

a a 

: drag coefficient 

C zero-lift drag coefficient 

DO 

C l : lift coefficient 

C : lift coefficient for maximum lift-to-drag ratio 

LR 

c : ratio of C /C 

L LR 

c : subscript for circularization or reorbit 

D : drag force 

d subsript for deorbit 
E : maximum of (L/D) 

m 

g gravitational acceleration 
g gravitational acceleration at surface level 

S 

H : altitude 

H : Hamiltonian 

h : subscript for Hohmann transfer 

i : subscript for idealized transfer 

J : performance index 

K induced drag factor 
L : lift force 

m : vehicle mass 

R : distance from vehicle center of gravity to Earth center 

R radius of the atmospheric boundary 
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R radius of the low Earth orbit 

C 

R : radius of the high Earth orbit 

d 

R radius of Earth 

E 

R distance from vehicle center of gravity to surface level 

S 

S : aerodynamic reference area 

s : subscript for surface level 

t : time 

V velocity 

v : normalized velocity 

£ inverse atmospheric scale height 

y : flight path angle 

S : normalized density 

A costate (Langrange) varible 

p : gravitational constant of Earth 

p : density 

r : normalized time 

AV : characteristic velocity 

Av : normalized characteristic velocity 

1. INTRODUCTION 

In space transportation systems, the concept of aeroassisted orbital 
transfer opens new mission opportunities, especially with regard to the 
initiation of a permanant space station [1]. The use of aeroassisted maneuvers 
to affect a transfer from high Earth orbit (HEO) to low Earth orbit (LEO) has 
been reccommended to provide high performance leverage to future space 
transportation systems. The space-based orbit transfer vehicle (OTV) is 
planned as a system for transporting payloads between LEO and other locations 
in space. The OTV, on its return journey from HEO, dissipates orbital energy 
through atmospheric drag to slow down to LEO velocity. In a synergetic 
maneuver for aeroassited orbital transfer vehicles (AOTV’s), the basic idea is 
to employ a hybrid combination of propulsive maneuvers in space and 
aerodynamic maneuvers in sensible atmosphere [2-4]. 


In this report, we first describe briefly the various types of coplanar 
transfers [5], Then we address the fuel-optimal control problem arising in 
coplanar orbital transfer employing aeroassist technology. The maneuver 
involves a transfer from HEO to LEO and at the same time minimization of the 
fuel consumption for achieving the desired orbit transfer. It is known that a 
change in velocity, also called the characteristic velocity, is a convenient 
parameter to measure the fuel consumption. A suitable performance index is the 
total characteristic velocity which is the sum of the characteristic 
velocities for deorbit and for reorbit (or circularization) [6]. Use of 
PontryaLgin minimum principle leads to a nonlinear, two-point, boundary value 
problem (TPBVP) in state and costate variables. This problem is solved by 
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using an efficient multiple shooting method [7,8]. Simulations are carried out 
for establishing a spectrum of entry conditions which are suitable for flying 
the spacecraft through the atmosphere. 

2. TYPES OF COPLAN AR ORBITAL TRANSFERS 

In transferring a vehicle from one orbit, say High Earth Orbit to another 
orbit, say Low Earth Orbit, in the same plane, we consider the following three 
types of orbital transfer. 

( i ) Hohmann transfer 

(ii) Idealized aeroassisted transfer 

(iii) Realistic aeroassisted transfer 

(i) Hohmann Transfer 

This is a coplanar transfer in which the vehicle is transferred from HEO 
to LEO in an all-proplusive manner (Figure 1). The first tangential impulse 
A V , called the deorbit impulse, at the radius R of HEO, is executed to 

d d 

transfer the vehicle from circular orbit to elliptic orbit with its perigee at 
R . The second tangential impulse AV , called the circularization or reorbit 

C C 

impulse, at the radius R of LEO, is impacted to transfer the elliptic orbit 

C 

to circular orbit. Thus, the total characteristic velocity AV h> a measure of 

fuel consumption in achieving the Hohmann transfer, is given by the sum of the 
deorbit impulse AV , and the cirularizat ion impulse AV . 

d c 

Using the principle of conservation of angular momenutum, the 
expressions for the deorbit and cirucularization impulses are given by [5] 



Using the normalized values, 

a = R /R ; a = R/R;Av = AV/[£7 r (3) 

d d A C C A A 

in (1) and (2), we get 

Av = TlTa. - <j2a /a (a +a ) (4) 

dh d c d d c 



so that the total characteristic velocity in the normalized form is- given by 
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(6) 


Av 

h 


Av ♦ Av 

dh ch 


Cii) Ideal Aeroassisted Transfer 

In an Idealized aeroassisted coplnar transfer, the vehicle grazes the 
atmospheric boundary, undergoes the necessary velocity reduction and skips 
back into another oriblt (Figure 2). Thus, the vehicle leaves HEO at R with 

d 

a tangential deorbit impulse AV and enters into an elliptic orbit with its 
perigee at R and flight path angle y = 0. When the vehicle is at perigee E, 

A C 

its lifting capability (in this case negative lift) is employed to affect 
flight path along the boundary of the atmosphere (i.e., along the circular 
orbit of radius R ). This graizing flight is continued along the atmospheric 

boundary to point F until sufficient velocity has been depleted by atmospheric 
drag such that upon reducing the lift to zero, the vehicle ascends with y = 0 

on an elliptic orbit to an apogee at R . Finally, at point C, a tangential 

C 

circularizing burn AV is imparted to achieve the desired LEO. In this 

C i 

transfer, the idealizations are (i) the atmospheric density at R^ is 

sufficient to generate enough drag to slow the vehicle in a reasonable amount 
of time, and (li) the vehicle has sufficient lifting capability to maintain 
fligh along the atmospheric boundany. 

The following relations for the characteristic velocities are obtained by 

using the principles of conservation of angular momenutum and the conservation 

of energy. The deorbit impulse AV given by 

dl 

AV = Jp/R- - 42(p/R )/(R /R ) (R /R +l) (7) 

dl d a ' d a d a 

and the circularization or reorbit impulse AV turns out to be 

ci 


AV = 4£^T - <l2(p/R )/(R /R )(R /R +T5 (8) 

cl c a'caca 

In terms of the normalized values (3), 

Av = 4l/a — 42/a (a +1) 

dl d d d 

Av = 4 1/a - 42/a (a +1) 

ci c c c 

and the total characteristic velocity for the idealized transfer in the 
normalized form is 


(9) 

( 10 ) 


Av = Av + Av 

l dl cl 


( 11 ) 
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(ill) Realistic Aeroassisted Transfer 


In a realistic, aeroassisted, coplanar transfer, the vehicle is 
transferred from HEO at R to LEO at R , by flying deep into the atmosphere to 

d c 

achieve the necessary velocity reduction (Figure 3). We start with a 
tangential propulsive burn, having a characterstic velocity AV for 

d 

deorbitting from the high Earth orbit (HEO) and entering into an elliptical 
transfer orbit. At point E the spacecraft enters the sensible atmosphere with 
an inclination of y and undergoes reduction in velocity due to atmospheric 
V e 

drag. At point F, the spacecraft leaves the atmosphere with flight path angle 
y . Once again, the transfer orbit is elliptical with the corresponding apogee 

at R . Finally, the maneuver ends with a cirularize or reorbit burn having a 

C 

characteristic velocity AV to maJce the vehicle enter into the low Earth 

C 

orbit. The desired circularization into LEO is achieved by the appropriate 
magnitude of AV . Thus, the maneuver consists of two impulses AV for deorbit, 

C d 

and AV for circularization. Upon applying the principles of conservation of 

C 

angular momentum and the conservation of energy, we obtain the relations for 
the deorbit and the circularization impulses as 


AV = 4p/R~ 1 

d d 


2( l-R /R )/fl-(R /R ) 2 /cos 2 y ) 

da 7 *- da' e' 


(12) 


AV = -Tp/fT 1 
c c 


2( l-R c /R a )/(l-(R c /R a ) 2 /cos 2 y f ) 


(13) 


In terms of the normalized values (3), 
Av = 4 1/a - J2(l-a )/a (l-a 2 /cos 2 y ) 

d d d ' d d e 

Av = 4 1/a - J2(l-a )/a (l-a 2 /cos 2 y ) 

c c c ' c c f 


(14) 


(15) 


Finally, the total characteristic velocity for the realistic aeroassisted 
transfer in the normalized form becomes 

Av = Av + Av (16) 

d c 

Let us note that if y = y =0, the relations (12) to (16) for the various 

e f 

characteristic velocities for the realistic transfer are the same ais those 
relations (7) to (11) for the corresponding characteristic velocities for the 
idealized transfer. Thus, the total characteristic velocity Av for the 
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realistic transfer is at least equal to that of Av^. In otherwords, the 
idealized transfer is the lower bound of the realistic transfer. 

The entry velocity v , and the exit velocity v f are also obtained in the 
normalized form as 


v = 2a (1-a )/(cos 2 y - a 2 ) 
e d d ' e d 


v = [2a ( 1-a )/(cos 2 y - a 2 ) 

t 4 c c ' f c 

(iv) Comparison of Orbital Transfers 

We now try to compare the three types of transfers, i.e., Hohmann 
transfer, ideal aeroasslsted transfer, and realistic aeroassisted transfer. 

For this, the following date is used. Radius of Earth, R e = 6,378 KM -, radius 

of atmosphere, R *= 6,498 KM; radius of the HEO, R j = 42,241 KM; radius of the 

a d 

LEO, R = 6,728 KM; gravitational constant of Earth, p = 3.986xl0 14 

C 

3 2 

meter /sec . Using this data, the deorbit and circularization velocities and 
the total characteristic velocities for each of the types of transfers is 
shown in Figure 4. The variation of entry velocity and the deorbit 

characteristic velocity AV with respect to entry flight path angle y are 

d e 

depicted in Figure 5. Also, the variation of exit velocity V and the reorbit 
characteristic velocity AV f are shown in Figure 6. The effect of altitude of 
LEO H , on the reorbit characteristic velocity AV is shown in Figure 7. As 

C C 

expected, the minimum AV at y increases with the altitude of LEO. The 

c c 

scenario at the entry and exit corridors shown in Figures 5 and 6 indicate 
that we need to fly the vehicle at near zero exit flight path angle for 
minimum fuel consumption. 

3. EQUATIONS OF MOTION 

Consider a vehicle with constant point mass, moving about a nonrotating 
spherical planet. The atmosphere surrounding the planet is assumed to be at 
rest, and the central gravitational field obeys the usual inverse square law. 
The equations of motion are given by, 

= Vsiny (19a) 

at 

^ = -AC V 2 exp(-H/3) - (p/R 2 )siny (19b) 

dt d 


(17) 

(18) 
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(19c) 


^ = ACVexp(-H/3) + [V/R - p/(R Z V)]cosy 

where A = Sp /2m , H = R-R , p = p exp(-H3) and 
8 8 8 

C = C + KC 2 for a drag polar. 

Using the normalized variables, 

t = t /1 r 2 /ji ; v = V/Vp/T 

and the dimensionless constants, 

h = H/H ; b = R M ; 6 = p/p = exp(-h/3H ); 

a a' a b a 

c = C l/ C l,- C l» = 


dh 

dr 


dv 

dr 


dy _ 

dr 


bvsiny 


- A b(l+c 2 )«v 2 - 


. , * ^ bvcosy 

V’ cdv * TFTThT 


b 2 slny 

(b-l+h) 2 

.2 

b cosy 
(b-l+h) 2 v 


where, A = C Sp H /2m; A = C Sp H /2m 

1 DO s a 2 LR 8 a 

4. OPTIMAL CONTROL 


( 20 ) 

(21a) 

(21b) 

(22a) 

(22b) 

(22c) 


For an optimal control problem with minimum fuel consumption, it is 
required to choose the performance index to minimize the total characteristic 
velocity, which is the sum of the initial characteristic velocity AV , the 

d 

deorbit impulse from HEO, and the final characteristic velocity AV , the 

C 

circularization impulse into LEO. Thus, the performance index is given by [6] 
J = AV = AV + AV (23) 

d c 

where AV = Vp/R - (R /R )V cos(-y ) 

d d a d e e 

AV = Vp/R - (R /R )V cosy 

c c a c f f 
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In the normalized form, the performance index becomes, 

J = Av = Av + Av (24) 

d c 

where , 

Av = vT7a - (v /a )cos(-y ) 

d ' d e d e 

Av = vT7a - (v /a )cos(y ) 

c • c t c f 


Alternatively, the expressions for AV and AV are given by (12) and 

d c 

(13). The minimization of J with respect to y and y f yields y = 0 and y = 0. 

Because of this, it is natural to postulate that the optimal trajectory 
behaves as y(t) = 0 for all values of time t between the entry time t and the 

e 

exit time t f . This constant flight path angle along the trajectory means that 

the altitude is constant at the value equal to altitude of the atmospheric 
interface. Thus, the optimal trajectory becomes an ideal (grazing) trajectory. 

We are interested in finding the minimization of the fuel with respect to 
the control C^. Using Pontryagin principle, we formulate Hamiltonian as 

H = A bvsiny + A Ia b(l+c 2 )Sv 2 - ~ b s * ny A 
h V t 1 (b-l+h) 2 J 


bvcosy 

(b-l+h) 


b 2 cosy * ) 
(b-l+h) 2 vJ 


where A^, A y , and A are the costates (Lagrange multipliers) corresponding to 
the states h, v, and y respectively. 

The unconstrained optimal control is obtained from 

|^ = 0 (26a) 

a c 

This leads us to 

c = E A /vA ; where, E = (L/D) (26b) 

m y' v o max 


Realistically, the control is bounded by the aerodynamic characteristics of 
the vehicle. Thus, for the contstrained control, 
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|C I s C or |c| s c 

1 L 1 Lmax < 1 ■ 


( 27 ) 


The costate (adjoint) variables A's are solved from 


dA 

h 

dr 


an 
a h’ 


dA 

v 

dx 


an 

av : 


dx dy 


( 28 ) 


Boundary Conditions 

The Initial and final boundary conditions are given for the normalized 
altitude h as 

h(x=0) = 1.0, h(x=x ) = 1.0 (29a) 

t 

and for the normalized velocity v, and the flight path angle y as 

(2-v 2 )a 2 - 2a - v 2 cos 2 y = 0 (29b) 

e d d e e 

(2-v 2 )a 2 - 2a - v 2 cos 2 y = 0 (29c) 

f c c f f 

The above relations result from the considerations of energy conservation 
and angular momentum conservation applied to the HEO-to-entry elliptic 
transfer orbit ard exit-to-LEO elliptic transfer orbit, respectively. The 
remaining boundary conditions are obtained from the transversal ity condi tons 
on the costates. Thus, the optimization procedure, requiring the solution of 
state equations (22) and the corresponding costate equations (28) along with 
the boudary conditons (29), leads us to the formation of a nonlinear, 
two-point, boundary value problem (TPBVP), which is solved by using a multiple 
shooting method. 

S. NUMERICAL DATA AND RESULTS 

A typical set of numerical values used for simulation purposes is given 
below [4], 

C =0.21; K = 1.67; m/S = 300 kg/m 2 

DO 

p =1.225 kg/m 3 ; p = 3.986xl0 14 m 3 /sec 2 

S 

0 = 1/6900 m ; R = 6378 KM 

E 

H = 120 KM; R = 12996 KM; R = 6558 KM 

ad c 

Using the above mentioned data, simulations are carried out. Time 
histories of altitude H, velocity V, ard flight path argle y, for total flight 
time of -540 seconds, are shown in Figure 8. Figure 9 shows the lift ard the 
lift-to-drag ratio ais a function of time from atmospheric entry to exit 
points. Those for the heating rate, dynamic pressure ard density are shown in 
Figure 10. The maximum positve L/D is used initially to recover from the 
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downward plunge. As the flight path angle becomes positve, the maximum 
negative L/D is used to level off the flight. Then a negative L/D is used to 
maintain flight at a small positve flight path angle in order to achieve the 
shallow exit. 

Figure 8(a) shows the time history of altitude. The spacecraft enters and 
exits the atmosphere at an altitude of 120 KM. The minimum altitude reached is 
55.56 KM. The velocity versus time is shown in Figure 8(b). The vehicle enters 
the atmosphere with a velocity of 9029 m/sec and leaves the atmosphere with a 
speed of 7795 m/sec, thus giving a velocity reduction of 1234 m/sec. The 
profile of flight path angle with time is shown in Figure 8(c). The spacecraft 
enters the atmosphere with an inclination of -5.665 degrees and exits with 
+0.927 degrees. The conrol history is shown is shown in Figure 9(a). The 
vehicle enters the atmosphere with maximum lift capability and switches to the 
minimum lift coefficient and then gradually increases during the ramaining 
flight. 

The minimum-fuel transfer requires a deorbit (impulse) characteristic 
velocity AV^ of 1034.29 m/sec and a reorbit characteristic velocity AV of 

73.25 m/sec, with a total characteristic velocity of 1107.54 m/sec. Let us 
compare this aeroassisted transfer with the Hohmann transfer, which is 
maneuvered entirely in outer space, and has a deorbit characteristic velocity 
AV of 1002.39 m/sec and the reorbit characteristic velocity AV of 1192.25 

dh ch 

m/sec, giving a total characteristic velocity for Hohmann transfer AV of 

h 

2194.64 m/sec. This shows that the saving due to coplanar, aeroassisted 
transfer over Hohmann transfer is 49.54 pecent. In the case of idealized 
transfer which follows a grazing trajectory along the atmospheric boundary, 
the deorbit characteristic velocity AV rfi is 1016.24 m/sec and reorbit 

characteristic velocity AV is 17.94 m/sec, thus giving a total 

C 1 

characteristic velocity of 1034.18 m/sec. The optimal transfer requires only 
6.63 percent more fuel than that of the idealized transfer. 

The heating rate Q , along the atmospheric trajectory, is computed for a 

r 

sphere of radius of 1 meter, according to the relation [3], 

Q = K p^V 3,08 (30) 

r r 

3 

where, p is the atmospheric density in kg/km , V is_the velocity in km/sec and 
K is the proportionality constant equal to 3.08x10 . Figure 10(a) shows peak 

r 

heating rate of 129.2 W/sq.cm. and the total integrated heating rate is found 
to be 15.536 KW-sec/cm . As shown in Figure 10(b), the peak dynamic pressure 
is 26.73 kN/sq.m. The variation of density in Figure 10(c) attains a maximum 
value of 0.3902 kg/m 3 . 

Entry Corridor 

A concept associated with atmospheric entry trajectories is that of the 
entry corridor. This is required due to the fact that a given vehicle cannot 
fly an acceptable atmospheric flight for orbitrary initial conditons at the 


11 



entry point. If the flight path angle y^ Is too steep, the vehicle will later 

suffer excessive aerodynamic and aerothermodynamic loadings even if the 
maximum lift is directed upward. This also may lead to "crash" condition. Or 
if the entry flight path angle is too shallow, the vehicle will exit the 
atmosphere again with an orbital velocity even if the maximum lift is directed 
downward. This leads to "escape" or uncontrolled skip-out condition. These 
boudaries of entry flight path angle are often taken to define the corridor of 
acceptable entry conditions. The responsibility of generating acceptable entry 
corridor lies largley with the last phase of the midcourse guidance system. 
Without a proper midcourse correction, it will be beyond the capability of the 
atmospheric flight control system to fly a safe trajectory. 

The entry corridor is the entry interface undershoot and overshoot and is 
usually specified by the entry flight path angle y^, as dictated by the entry 

dynamics. In the present case of fuel-optimal, coplanar, orbital transfer, 
four simulations are carried out as shown in Table I below. Figure 11 shows 
the entry corridor with flight path angle ranging from -7.24 deg -5.485 
degrees. 

Table I: 


Simulation 

No 

r 

e 

deg 

V 

e 

m/sec 

C L 

flight time 
seconds 

1 

-7.240 

9020 

0.5381 

510 

2 

-6.412 

9025 

0.5786 

600 

3 

-5.665 

9029 

0.6299 

540 

4 

-5. 485 

9030 

0.6502 

600 


Approximate Solutions 

As a first cut in getting approximate solutions, we try to obtain the 
optimal solutions using the fact that the flight path angle y is small. Using 
this approximation, the simulation has been repeated. Interestingly, it is 
found that we get almost the identical solutions with much less computation 
time. For example, in this case, the deorbit characteristic velocity AV is 

d 

equal to 1034.29 m/sec, and the reorbit characteristic velocity AV is found 

C 

to be 73.26 m/sec, giving a total minimum characteristic velocity as 1107.55 
m/sec. The total integrated heating rate is 15.556 KW-sec/cm . The execution 
time on the computer is now only 59.2 percent of the execution time for the 
optimal solution without the approximation. 

6. MULTIPLE SHOOTING METHOD 

The determination of optimal control (26) requires the solution of a 
sixth order, nonlinear TPBVP consisting of state equations (22) and costate 
equations (28) and the associated boundary conditions (29). This can only be 
done by numerical methods. The multiple shooting method is one of the 
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powerful methods for solving nonlinear TPBVP’s. The corresponding OPTSOL code 
was developed by DFVLR establishment at Oberpfaffenhofen, West Germany [7,8]. 

In solving any boundary value problem with the given initial and final 
conditions, we asume additional initial data and integrate forward so that the 
solution satisfies the given final condition as well. This is also called a 
simple shooting method. Here, the convergence of the solution is highly 
sensitive to the assumed intial data. It is found that the error due to 
inaccurate initial data can be made arbitrarily small by performing the 
integration over sufficiently smaller subdivided panels within the given 
interval and thereby leading to the multiple shooting method. Thus, the 
multiple shooting method is a simultaneous application of the simple shooting 
method at several points within the interval of integration. Here, the 
trajectory may be restarted at intermediate points using new guesses. Jacobian 
matrices are formed for each segment. The resulting iteration scheme, based on 
reducing all discontinuties at internal grid points to zero, leads to a system 
of linear algebraic equations. 

Figure 12 shows the successive approximations of the altitude H, during 
the course of 0, 5, and 14 iterations. For the sake of clarity only 4 out of 
20 intervals are shown. The initial guessed value for the altitude is 120 KM 
at every interval. It can be seen how the intitially large jumps at the 
subdivision points of the multiple shooting method are "flattened out" with 
the increase of iterations. 

7. CONCLUSIONS 

A breif treatment of the three types of orbital transfers, Hohmann 
transfer. Idealized transfer, and realistic transfer, established the fact 
that the realistic transfer is the upper bound of the idealized (grazing) 
transfer. 

In this report, we have addressed the mlmimization of fuel consumption 
during the atmospheric portion of an aeroassisted, coplanar, orbital transfer 
vehicle. The resulting two-point, boundary value problem was solved by using 
an efficient multiple shooting algorithm. Simulations have been carried out to 
establish a corridor of entry conditions suitable for flying the vehicle 
through the atmosphere. The strategy for the atmospheric portion of the 
minimum-fuel transfer is to fly at the maximum L/D initially in order to 
recover from the downward plunge, and then to fly at a negative L/D to level 
off the flight such that the vehicle skips out of the atmosphere with a flight 
path angle near zero degrees. 
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Fig.l Hohmann transfer 
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Fig. 2 Ideal aeroassisted transfer 
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Fig. 4 Comparison of orbital transfers 
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Fig. 6 Relations at the exit corridor 
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Fig. 7 Effect of LEO on reorbit characteristic velocity 
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Fig. 8(c) Time history of flight path angle 
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Fig. 10(a) Time history of heating rate 
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Fig. 10(c) Time history of density 
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Fig. 12 Successive approximations for altitude 


FUEL-OPTIMAL TRAJECTORIES FOR AERO ASS IS TED 
NONCOPLAN AR ORBITAL TRANSFER PROBLEM 


Dr. D. S. Naidu 
ODU Research Foundation 
Norfolk, VA, 23508 


ABSTRACT: The fuel-optimal control problem arising in noncoplanar orbital 
transfer employing aeroasslst technology is addressed. The maneuver involves 
the transfer from high Earth orbit to low Earth orbit with plane change. A 
performance index is chosen to minimize the total fuel consumption for the 
transfer. Use of Pontryagin minimum principle leads us to a nonlinear, 
two-point, boundary value problem, which is solved by using a multiple 
shooting method. 


NOMENCLATURE 

A constant = Sp /2m 

B 7 

A : constant = C Sp H /2m 

1 dO s bf 

A : constant = C Sp H /2m 

2 LR s tJ 

b : ratio of R Ai 

a' a 

C d : drag coefficient 

C zero-lift drag coefficient 

DO 

C : lift coefficient 

L 

C : lift coefficient for maximum lift-to-drag ratio 

LR 

c : ratio of C /C 

L LR 

c : subscript for circularization or reorbit 

D : drag force 

d : subsript for deorbit 

E : maximum of (L/D) 

m 

g : gravitational acceleration 

g : gravitational acceleration at surface level 

S 

H : altitude 

3f : Hamiltonian 

J : performance index 

K : induced drag factor 

L : lift force 

m : vehicle mass 

R : distance from vehicle center of gravity to Earth center 

: radius of the atmospheric boundary 

R : radius of the low Earth orbit 

C 

R^ : radius of the high Earth orbit 

R : radius of Earth 

E 

R : distance from vehicle center of gravity to surface level 

S 

S : aerodynamic reference area 
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s : subscript for surface level 

t : time 

V : velocity 

v : normalized velocity 

0 : inverse atmospheric scale height 

y flight path angle 

0 : heading angle 

cr bank angle 

0 : down range angle or longitude 

<f> cross range angle or latitude 

5 : normalized density 

A : costate (Langrange) varible 

H : gravitational constant of Earth 

p : dens i ty 

x : normalized time 

AV : characteristic velocity 

Av : normalized characteristic velocity 

1. INTRODUCTION 

In space transportation system, the concept of aeroassisted orbital 
transfer opens new mission opportunities, especially with regard to the 
initiation of a permanant space station [1]. The use of aeroassisted maneuvers 
to affect a transfer from high Earth orbit (HEO) to low Earth orbit (LEO) has 
been reccommended to provide high performance leverage to future space 
transportation systems. The space-based, orbit transfer vehicle (OTV) is 
planned as a system for transporting payloads between LEO and other locations 
in space. The OTV, on its return journey from HEO, dissipates orbital energy 
through atmospheric drag to slow down to LEO velocity. In a synergetic 
maneuver for aeroassited, orbital transfer vehicles (AOTV’s), the basic idea 
is to employ a hybrid combination of propulsive maneuvers in space and 
aerodynamic maneuvers in sensible atmosphere. Within the atmosphere, the 
trajectory control is achieved by means of lift and bank angle modulations 
Hence, this type of flight with a combination of propulsive and nonpropulsive 
maneuvers, is also called synergetic maneuver or space flight[2-7]. 

In this report, we address the fuel-optimal control problem arising in 
noncoplanar orbital transfer employing aeroassist technology. The maneuver 
involves the transfer from HEO to LEO with a presribed plane change and at the 
same time minimization of the fuel consumption. It is known that the change in 
velocity, also called the characteristic velocity, is a convenient parameter 
to measure the fuel consumption. For minimum-fuel maneuver, the objective is 
then to minimize the total characteristic velocity for deorbit, boost, and 
reorbit (or circularization). Use of Pontryagin minimum principle leads us to 
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a nonlinear, two-point, boundary value problem (TPBVP), which is solved by 
using a multiple shooting method [10-12]. 

2. BASIC EQUATIONS 

The basic equations to be used are those for deorbit, aeroassist (or 
atmospheric flight), boost and reorbit (or circularization). 

(i) Deorbit 

Initially, the spacecraft is in cirualar orbit of radius R , well outside 

d 

the Earth’s atmosphere, moving with a circualar velocity V = Vfi/R . Deorbit 

d d 

is accomplished by means of an impulse AV , to transfer the vehicle from a 

d 

circular orbit to elliptic orbit with perigee low enough to intersect the 
dense part of the atmosphere. Since the elliptic velocity at D is less than 
the circular velocity at D, the impulse AV is executed so as to oppose the 

d 

circular velocity V . In otherwards, at point D, the velocity required to put 

d 

the vehicle Into elliptic orbit is less than the velocity required to maintain 
it in circular orbit. The deorbit impulse AV causes the vehicle to enter the 

d 

atmosphere at radius R with a velocity V and flight path angle y . It is 

A - C © 

known that the optimal -energy loss maneuver from the circular orbit is simply 
the Hohmann transfer and the impulse is parallel and opposite to the 
instantaneous velocity vector. 

Using the principle of conservation of energy and angular momentum at the 
deorbit point D, and the atmospheric entry point E, we get [8], 

V 2 /2 “ p/R = (V -AV ) 2 /2 - p/R (1) 

R V cos(-y ) = R (V -AV ) (2) 

a e e add 

from which solving for AV , we get 

d 

AV = 4p/R" - 2p( 1/R -1/R )/r(R /R ) 2 /cos\ -ll (3) 

d d H a d ' L d a ' e J 


It is easily seen that the minimum value of the deorbit impulse AV 

dm 

obtained at y =0, corresponds to an ideal transfer with the space vehicle 

grazing the atmospheric boundary. To ensure proper atmospheric entry, deorbit 
impulse AV must be higher than the mimimura deorbit impulse AV which is 

d dm 

given by 
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( 4 ) 


AV = Jp/R - - |2fi(l/R -1/R )/f(R /R ) 2 -l1 

dm d *i a d ' d a J 

(ii) Aeroassist 

During the aeroassist (or atmospheric flight), the vehicle needs to be 
controlled by lift and bank angle to achieve the necessary velocity reduction 
(due to atmopheric drag) and the plane change. 

Consider a vehicle with constant point mass m, moving about a nonrotating 
spherical planet. The atmosphere surrounding the planet is assumed to be at 
rest, and the central gravitational field obeys the usual inverse square law. 
The equations of motion for the vehicle are given by (Figure 1), 


dH ,, , 

at = Vslny 

(5a) 

~ = - AC V 2 exp(-H/3) - ( p/R 2 ) s 1 ny 
dt d 

(5b) 

^ = AC L Vcos<rexp(-H/3) + [V/R - p/(R 2 V)]cos y 

(5c) 

d6 

= Vcosycosi/»/(Rcos0) 

(5d) 

^ = (V/R)cosysin^ 
at 

(5e) 

^ = AC L Vsin<rexp(-H0)/cosy - (V/RJcosycos^tan# 

(5f ) 


where A = Sp /2m , H = R-R , p = p exp(-H0) and 

S S 8 

C = C +■ KC 2 for a drag polar. 

D DO L “OK 

Using the normalized variables, 

t = t/!5> * v = V/Vp/R (6) 

and the dimensionless constants, 

h = H/H ; b = R/H; 8 = p/p = exp(-h£H ) 

a a' a s a 

c = C /C ; C = {C /K 

V LR LR * DO' 


(7a) 

(7b) 
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in (5), we get the normalized form as 


dh 

dx 


bvsiny 


(8a) 


dv 

dx 


2 

- Ab(l + cV - — S -ig 

1 (b-l+h) 2 


dy 

dx 


A bcSvcosc 
2 


bvcosy 

(b-l+h) 


.2 

b cosy 
(b-l+h) 2 v 


d9 _ bvcosycos^ 
dx ~ (b-l+h)cos0 


(8b) 

(8c) 

(8d) 


d <p _ bvcosysin^ 
dx = (b-l+h) 

, . A b5Avsin<r , .. 

dip _ _2 bvcosycos0taji0 

dx ~ cosy (b-l+h) 


(8e) 


(8f ) 


where, A = C Sp H /2m; A = C Sp H /2m 

1 DO s a 7 2 LR K s a 7 


The plane change or orbit inclination, i, is related with cross range 
angle <p and heading angle \p as [9] 

cosi = cos<f>cosil> t^tst^. (9) 

The orbit inclination changes through out the atmosphere segment and must have 
the required value at exit. For small values of cross range angle <p, i is 
given by the heading angle ip Itself. 

(iii) Boost and Reorbit 

During the atmospheric flight, the vehicle undergoes the plane change due 
to modulation in lift and bank angle. Because of the loss of energy during a 
turn, a second Impulse is required to boost the vehicle back to orbital 
altitude. 

The vehicle exits the atmosphere at point F, with a velocity V and 

flight path angle y . The additional impulse AV , required at the exit point F 

f b 

for boosting into an elliptic orbit with apogee radius R and the reorbit (or 

C 

circularization) impulse AV required to insert the vehicle into a circular 

C 

orbit, are obtained by using the principle of conservation of energy and 
angular momentum at the exit point F, and the reorbit or circularization point 
C. Thus, we have, 
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(V +AV ) 2 /2 - p/R = (V -AV ) Z /2 - p/R 

f b ' a c c ' c 


( 10 ) 


(V +AV )R cosy » R (V -AV ) 

f b a f c c c 

Solving for AV and AV from the above equations (10) and (11), 

b c 


(id 


AV = 2p( 1/R -1/R )/ri-(R /R ) 2 cos 2 y 1 - V 

b i a c ' «■ ac f 


( 12 ) 


AV = - 2p( 1/R -1/R )/[(R /R ) 2 /cos 2 y -l] 

C C *1 A C “ C & » * 


(13) 


3. OPTIMAL CONTROL 

It is known that the change in speed, AV, also called the characteristic 
velocity, is a convenient parameter to measure the fuel consumption. For 
minimum-fuel maneuver, the objective is then to minimize the total 
character i stic velocity. A convenient performance index is the sum of the 
characteristic velocities for deorbit, boost, and reorbit. Thus, 

J = AV + AV + AV (14) 

d b c 

Where, A V , AV , and AV are the deorbit, boost, and reorbit characteristic 

d b c 

velocities respectively, and are related as 


AV = Vp/R - (R /R ) V cos(-y ) 

d d a d e e 


(15) 


AV = Vp/lT - (R /R ) ( V +AV )cosy (16) 

c c a c f b f 


Alternatively, AV , AV and AV are also given by (3), (12), and (13) 

d b c 


respectively. 

In the normalized form, the performance index becomes. 


J = Av = Av 

+ Av + Av 

(17) 

d 

where , 

b c 

Av = Vl /a • 

d ' d 

- (v /a )cos(-y ) 

e d e 

(18) 

Av = /l7a ■ 

c 7 c 

- [(v +Av ) /a Icosy 

L f b 9 c J f 

(19) 
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( 20 ) 


a d = W % = W’ AV = AV /^ a ' 

Let us note that for a given circular orbit, the Impulses AV and AV are 

b c 

completely determined by the state variables V f and y at the exit conditions 
of the atmospheric portion of the trajectory. The velocity V^ and the flight 
path angle y^ at the entry point are dependent only on the magnitude of the 
deorbit impulse AV . Therefore, the optimal control problem needs to consider 

d 

the segment of the trajectory within the atmosphere. 

The first step in the optimization procedure using Pontryagin principle 
is to formulate Hamiltonian els [9] 


H = 


A bvsiny 
h 


+ A 


-t 


A b(l+c 2 )«v 2 


b 2 siny 1 
(b-l+h) 2 J 


+ A 


r{‘ 


A bc5vcos<r + 
y 2 


bvcosy b 2 cosy "I 

rerTTKT_ (b-l+h) 2 vJ 


- fbvcosys i ni/Tl 

+ X <f>\ (b-l+h) j 


+ A 


& 


A 2 b6Avsin(r bvcosycos^tan^l 


cosy 


(b-l+h) J 


( 21 ) 


where A’s are the costates corresponding to the states. The down range angle 0 
does not enter the right hand side of the equations of motion (5) and hence 
need not be considered for the optimization process. 

The optimal control equations for lift and bank angle Eire given by 


an 

dc 


- 0 ; 



( 22 ) 


leading to 

c = C w/C 2vA ; tan<r = A, A cosy 

LR 'DO v if/ y 


(23) 


where 


W * j*y + CV COSy 3 ! 


(24) 


✓ 
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Realistically, the control is bounded by the aerodynamic characteristics of 
the vehicle. Thus, for the constrained control, 


|C I s C or |c| s c 

1 L ' Lmax * * max 

The costate (adjoint) equations are given by 


dA 

h 

an 

dA 

V 

dx 

~ ah’ 

dx 


dK 

dA. 

_ Jl 

dx 

" d <t> ; 

dx 


8H 

3v ; 


dH 

d\ft 



8H 

87 


(25) 


(26a) 


(26b) 


Boundary Conditions 

The initial and final boundary conditions are given for the normalized 
altitude h as 

h(x=0) = 1.0, h(x=X f ) = 1.0 (27a) 

and for the normalized velocity v, and the flight path angle 7 as 

(2-v Z )a 2 - 2a - v 2 cos 2 * =0 (27b) 

e d d e e 

T2-(v +Av ) 2 la 2 - 2a - (v +Av ) 2 cos 2 r = 0 (27c) 

•-fb J ccfb f 

The above relations are obtained using the principle of energy conservation 
and angular momentum conservation to the HEO-to-entry elliptic transfer orbit 
and the exit-to-LEO elliptic transfer orbit, respectively. The remaining 
boundary conditions are obtained from the transversal ity conditions on the 
costates. Thus, the optimization procedure, requiring the solution of the 
state equations (8) aind the costate equations (26) along with the boundary 
conditions (27) leads us to the formulation of a nonlinear TPBVP, which is 
solved by using a multiple shooting method [10-12]. 

4. NUMERICAL DATA AND RESULTS 

A typical set of numerical values used for simulation purposes is given 
below [3-6], 

C =0.1; K = 1.111; m/S = 300 kg/m 2 

DO 

p =1.225 kg/m 3 ; p = 3.986xl0 14 m 3 /sec 

S 

B = 1/6900 m" 1 ; R £ = 6378 KM 
H = 120 KM; R = 12996 KM; R = 6558 KM 

a d c 
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Using the above mentioned data, simulations are carried out. The optimal 
solution has the following entry and exit status. 

Entry status: H = 120 KM; V = 9034.74 m/sec 

e e 

7 * = -4.435 degrees; 0^ = 0; 0^ = 0 

Exit status: H = 120 KM; V = 7009.29 m/sec 

f f 

7 = -0.6217 deg; 0 f = 7.237 deg 

0 f = 18.467 deg; total flight time = 500 sec 

Characteristic velocities: 

Deorbit characteristic velocity, AV = 1034.29 m/sec 

d 

Boost characteristic velocity, AV = 816.00 m/sec 

b 

Reorbit characteristic velocity, AV = 42.96 m/sec 

C 

Total characteristic velocity AV = 1893.25 m/sec 

Time histories of altitude H, velocity V, and flight path angle 7 , for 
total flight time of 540 seconds, are shown in Figure 2. Figure 3 shows lift 
C^, bauik angle <r, and lift-to-drag ratio L/D as a function of time from 

atmospheric entry to exit points. Those for the heating rate, dynamic pressure 
and density are shown in Figure 4. 

Figure 2(a) shows the time history of altitude. The spacecraft enters and 
exits the atmosphere at an altitude of 120 KM. The minimum altitude reached is 
46 KM. The velocity versus time is shown in Figure 2(b). The vehicle enters 
the atmosphere with a velocity of 9034.74 m/sec and leaves the atmosphere with 
a speed of 7009.29 m/sec, thus giving a velocity reduction of 2025.45 m/sec. 
The profile of flight path angle with time is shown in Figure 2(c). The 
spacecraft enters the atmosphere with an inclination of -4.435 degrees and 
exits with -0.6217 degrees. The conrol history is shown is shown in Figure 
3(a). The vehicle enters the atmosphere with maximum lift capability and 
decreases slowly during the ramaining flight. Figure 3(b) shows the variation 
of bank angle during the atmospheric flight. Initially the vehicle enters the 
atmosphere with a bank angle of 154.75 degrees to pull the vehicle into the 
atmosphere but slowly drops to about 63.3 degrees auid maitain at that value 
for the remainder of the flight. 

The minimum-fuel transfer requires a total characteristic velocity of 
1893.25 m/sec. It is of interest to compare this optimal aeroassisted transfer 
with the Hohmann transfer, which is maneuvered entirely in outer space, and 
has a deorbit characteristic velocity AV of 1002.39 m/sec and the reorbit 

dh 

characteristic velocity AV of 1192.25 m/sec, giving a total characteristic 

ch 

velocity for Hohmann transfer AV of 2194.64 m/sec. This shows that still the 

h 

noncoplanar transfer with a plane change of 18.467 degrees is only 84 percent 
of the Hohmann (coplanar) transfer. 
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The heating rate Q , along the atmospheric trajectory, is computed for a 

r 

sphere of radius of 1 meter, according to the relation [6], 

Q = K p^V’ 08 (28) 

r r 

3 

where, p is the atmospheric density in kg/km , V is the velocity in km/sec and 
K is the proportionality constant equal to 3. 08x10 . Figure 4(a) shows the 

r 

peak heating rate of 251.16 W/sq. cm. and the total integrated heating rate is 
found to be 21.0435 KW-sec/cm . As shown in Figure 4(b), the peak dynamic 
pressure is 103 KN/sq. m. The variation of density in Figure 4(c) ottained a 
peak value of 1.5367 kg/m 3 . 

Approximate Solutions 

As a first cut in getting approximate solutions, we try to obtain the 
optimal solutions using the fact that the flight path angle y is small. Using 
this approximation, the simulation has been repeated. Interestingly, it is 
found that we get almost the identical solutions with much less computation 
time. For example, in this case, the deorbit characteristic velocity AV is 

d 

equal to 1034.29 m/sec, boost characteristic velocity AV is 816 m/sec and the 

b 

reorbit characteristic velocity AV is found to be 43.01 m/sec, giving a total 

C 

minimum characteristic velocity as 1893.3 m/sec. The total integrated heating 
rate is 21.0435 KW-sec/cm . The execution time on the computer is now only 41 
percent of the execution time for the optimal solution without the 
approx i mat ion. 

6. MULTIPLE SHOOTING METHOD 

The determination of optimal control (23) requires the solution of a 
tenth order, nonlinear TPBVP consisting of state equations (8) and costate 
equations (26) and the associated boundary conditions (27). This can only be 
done by numerical methods. The multiple shooting method is one of the 
powerful methods for solving nonlinear TPBVP’ s. The corresponding 0PTS0L code 
was developed by DFVLR establishment at Oberpfaffenhofen, West Germany 
[ 10 - 12 ]. 

In solving any boundary value problem with the given initial and final 
conditions, we asume additional Initial data and integrate forward so that the 
solution satisfies the given final condition as well. This is also called a 
simple shooting method. Here, the convergence of the solution is highly 
sensitive to the assumed intial data. It is found that the error due to 
inaccurate initial data can be made arbitrarily small by performing the 
integration over sufficiently smaller subdivided panels within the given 
interval and thereby leading to the multiple shooting method. Thus, the 
multiple shooting method Is a simultaneous application of the simple shooting 
method at several points within the interval of integration. Here, the 
trajectory may be restarted at intermediate points using new guesses. Jacobian 
matrices are formed for each segment. The resulting iteration scheme, based on 
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reducing all discontinuties at internal grid points to zero, leads to a system 
of linear algebraic equations. 

Figure 5 shows the successive approximations of the altitude H, during 
the course of 0, 5, and 14 iterations. For the sake of clarity only 4 out of 
20 intervals are shown. The initial guessed value for the altitude is 120 KM 
at every interval. It can be seen how the intitlally large jumps at the 
subdivision points of the multiple shooting method are "flattened out" with 
the increase of iterations. 

7. CONCLUSIONS 

In this report, we have addressed the mimimization of fuel consumption 
during the atmospheric portion of an aeroassisted, noncoplanar, orbital 
transfer problem. The resulting two-point, boundary value problem was solved 
by using an efficient multiple shooting algorithm. The strategy for the 
atmospheric portion of the minimum-fuel transfer is to start initially with 
the maximum postive lift in order to recover from the downward plunge, and 
then to fly with a gradually decreasing lift such that the vehicle skips out 
of the atmosphere with a flight path angle near zero degrees. Also, initially, 
a bank angle greater than 90 degrees is used to pull the vehicle into the 
atmosphere, but is later reduced to produce the skip. 
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Fig. 2(a) Time history of altitude 
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Fig. 2 (b) Time history of velocity 
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Fig. 2(c) Time history of flight path angle 
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Fig. 2(d) Time history of cross range 
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Fig. 2(e) Time history of heading angle 
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Fig. 3(b) Time histo^ of bank angle 
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Fig. 4(b) Time history of dynamic pressure 
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Fig. 5 Successive approximations for altitude 


